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Abstract 

The order-disorder duality structure is exploited in order to ob- 
tain a quantum description of anyons and vortices in: a) the Maxwell 
theory; b) the Abelian Higgs Model; c) the Maxwell-Chern-Simons 
theory; d) the Maxwell-Chern-Simons-Higgs theory. A careful con- 
struction of a charge bearing order operator(a) and a magnetic flux 
bearing disorder operator (vortex operator) (/i) is performed, paying 
attention to the necessary requirements for locality. An anyon opera- 
tor is obtained as the product (p = cru. A detailed and comprehensive 
study of the euclidean correlation functions of a, \l and ip is carried 
on in the four theories above. The exact correlation functions are 
obtained in cases a and c. The large distance behavior of them is 
obtained in cases b and d. The study of these correlation functions 
allows one to draw conclusions about the condensation of charge and 
magnetic flux, establishing thereby an analogy with the Ising model. 
The mass of vortex and anyon excitations is explicitly obtained wher- 
ever these excitations are present in the spectrum. The independence 
between the mechanisms of mass generation for the vortices and for 
the vector field is clearly exposed. 
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1) Introduction 



The physics of 2 + 1 dimensional world has been the object of intense theoretical 
investigation in the last few years. One of the features which has aroused a lot of 
interest is the possibility of existence of states with generalized, continuous statistics. 
Chern-Simons field theories [1] have been playing a central role in this framework. 
An important reason for this is that the presence of a Chern-Simons term in the 
lagrangian of a vector field induces a change in the statistics of the charged particles 
eventually coupled to it [2]. This statistical transmutation effectively occurs because 
the coupling of a charged particle to a Chern-Simons vector field imparts to this 
particle a point magnetic flux [2] and, as was demonstrated in full generality [3], the 
bound state of a point charge and a point magnetic flux carries arbitrary statistics 
proportional to the product: charge x magnetic flux. A particular case of this fact 
was first observed in a specific model in [4]. 

The mechanism of statistical transmutation via the Chern-Simons lagrangian is a 
key ingredient of field theories proposed to describe very interesting condensed matter 
systems such as the two dimensional electron gas undergoing the Quantum Hall effect 
[5]. It was also invoked in the theory of superconductivity [6] involving particles with 
generalized statistics or anyons as they become known. 

In view of their vast potential physical interest, a large amount of study of anyon 
properties was undergone recently [7]. In order to obtain a full quantum field theoretic 
description of anyons one should be able to evaluate their correlation functions, mass 
spectrum and scattering amplitudes. The purpose of this work is to accomplish some 
of these goals in the framework of four different quantum field theories: a) the Maxwell 
theory; b) the Abelian Higgs Model (Maxwell-Higgs theory); c) the Maxwell-Chern- 
Simons theory; d) the Maxwell-Chern-Simons-Higgs theory. 

A lot of work has been devoted to the study of generalized statistics in 1 + 1 D 
[8,4,9]. Actually the investigation performed in 1 + ID provided a fundamental clue 
for understanding the quantum description of anyons, either in 1 + ID and 2 + ID. 
This consists in the realization that there is a basic algebraic structure underlying 
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the existence of generalized statistics excitations, namely, the order- disorder duality 
[10]. The idea is to introduce order (a) and disorder (//) operators satisfying a certain 
algebra (dual algebra) in such a way that the composite operator <p — a/i made out 
of them is shown to be in general an anyon operator [4]. One then generalizes the 
method of Kadanoff and Ceva [10] for continuum field theory [4,11] in order to obtain 
the euclidean correlation functions of a, ji and <p. A review of this method is in ref. 
[12]. Very interesting and important related works are found in [13]. 

A few years ago, the order- disorder duality structure was generalized for 2 + 1- 
dimensional continuum field theory [14]. A magnetic flux bearing disorder operator \x 
was introduced and its correlation functions obtained as euclidean functional integrals. 
It was also shown that the order operator a dual to ji should now be a charge bearing 
operator. As a consequence, the product cr/i is naturally seen to be an anyon operator. 
This description of anyons goes beyond the Chern-Simons mechanism because one 
can study generalized statistics excitations even in theories without a Chern-Simons 
term. Very interesting related works were done previously on the the lattice [15] and 
subsequently in the continuum [16]. 

An extremely interesting and appealing feature of this approach emerges when one 
realizes that the magnetic flux is the topological charge in 2 + ID. Hence the disorder 
operator /i or the anyon itself, must be closely related to topological charge bearing 
excitations in 2 + ID, namely, topological solitons or vortices. This relationship was 
investigated and demonstrated in [14,16] where the disorder variable ji was shown to 
be a quantum vortex creation operator. 

The above ideas were recently applied in the complete bosonization of the Dirac 
fermion field in 2 + ID [18]. In this case the a and /i dual operators were expressed 
in terms of a bosonic vector field with dynamics governed by a nonlocal version 
of the Maxwell action plus a Chern-Simons term. 

It follows from the work done in [14] that some stringent requirements must be met 
in order for the vortex operator to be local. Also, a lot of care must be exercised in 
order to obtain the charge bearing local order operator a dual to it. We will carefully 
analyze the construction of local order, disorder and anyon operators in the various 
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theories under consideration, according to the requirements set forth in [14]. 

After expressing the fields a, /i and (p in terms of a vector field W^, we make a 
detailed study of the euclidean correlation functions of these operators and evaluate 
the commutation rules among themselves as well as between them and the charge 
and magnetic flux operators in the four theories mentioned above. Exact results for 
all correlation functions are obtained in the case of the Maxwell and Maxwell-Chern- 
Simons theories. The long distance behavior of the correlation functions is obtained 
in the case of the Abelian Higgs Model and Maxwell-Chern-Simons-Higgs theories, 
both in the symmetric and broken phases. 

We investigate the conditions which allow the presence of genuine vortex on anyon 
excitations in the spectrum and in these cases, evaluate their masses, for all of the 
above mentioned theories. 

As in 1 + 1D [4], we observe that a nontrivial commutation relation between opera- 
tors, manifests itself in the multivaluedness of the corresponding euclidean correlation 
functions. Analytic continuation back to Minkowski space, from each sheet, will lead 
to the various orderings of operators. We are thereby able to retrieve the commuta- 
tion rules directly from the study of euclidean correlation functions. We show that 
the general expression for the statistics of a field bearing a charge Q and a magnetic 
flux $ is S=^ , where Q = Q — 0<E> with 9 being the coefficient of the Chern-Simons 
term eventually present in the theory. 

We can envisage some interesting applications of our formalism in field theories 
related to condensed matter systems. Among them we could mention the Chern- 
Simons-Landau-Ginsburg theory of the Quantum Hall Effect [19] and the theories 
describing the Hall-Superconducting phase transition in anyon systems [20]. We are 
presently investigating such applications. 

The material is organized in the next four chapters in such a way that each chapter 
corresponds to one of the theories mentioned above. Concluding remarks and the main 
conclusions are presented in chapter 6. Six Appendixes are included to demonstrate 
useful results. 
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2) The Maxwell Theory 
2.1) Order and Disorder (Vortex) Operators 

2.1.1) The Disorder (Vortex) Operator 

Let us start by considering the Maxwell theory, described by 

C = -\w^W^ (2.1) 

which leads to Maxwell equation 

d v W v ' i = (2.2) 

The identically conserved topological current is given by 

= ^daWf, (2.3) 

We can also introduce the two indexes dual topological current 

J^ = e^ a J a or j^ = W^ (2.4) 

is not identically conserved; rather, its conservation is implied by the field equation 
(2.2). 

A topological charge (magnetic flux) bearing operator (vortex operator) was first 
introduced as a disorder variable in continuum 2 + 1 D field theory in [14]. A fun- 
damental piece in the construction of the vortex operator in [14] and which was also 
used in [18] is the singular external field 

A,(z; x- C- T x ) = f d%arg(i- x)5 3 (z - f) (2.5) 

JT X {C) 

In this expression, arg(z) is the angle the vector z makes with the z 1 axis (z = (z°, zj). 
T X (C) is the portion of the R 2 plane at £° = x°, external to the curve C. C is the curve 
depicted in Fig. 1 which contains the arc of a circumference of radius p, centered on 
(x°,x) and two straight lines along the cut of arg(z-x) separated an angle 25 apart. 
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The vortex operator (disorder operator) /i is one of the two basic dual operators 
whose properties we are going to investigate in this work. It is obtained by coupling 
the external field to the dual current J^ u = W^ u in the following way [14] 

H(x; C) = exp{~ J d 3 zW^A, u } (2.6) 

In this expression, b is a free real parameter with dimension of (mass)~^ and = 
d^Ay — djjA^. In [14] it is shown that the correlation functions of \i are completely 
independent of the surface T x appearing in the definition of A^, provided we introduce 
a surface renormalization counterterm depending on A^ and having the same form as 
the kinetic term of the field, namely 

Cr = ( 2 - 7 ) 

In the presence of this counterterm, it is clear that the integrand in the functional 
integral defining /i correlation functions depends on the external field A^ through the 
conbination + A^.( See Eq.(2.26)) This is a fundamental requirement for surface 
invariance because making a change of variable of the type — > W = + d^cu, 
which of course leaves the measure invariant we have W^ + A^ — > W + A^ — d^cu . One 
can demonstrate [14] that for an appropriate choice of u, we get A^(z; x; C; T x )—d^u = 
^(zjzjC;^), where 

A„(z;x;C;f x ) = f d%arg{£ - x)5 3 (z - £) - Q(V(f x ))d^arg(z - x) (2.8) 

Here T x is an arbitrary surface bounded by C and V(T X ) is the volume enclosed by 
T X (C) and T X (C). Q(V(T X )) is the three-dimensional heaviside function with support 
on V(f x ). (See Fig. 2 for f x (C) and V{f x )). It is not difficult to realize [14] that the 
uj which leads from (2.5) to (2.8) is 

uj = Q{V{f x ))arg{z-x) (2.9) 

A^z; x; C; T x ), eq. (2.8), is the most general form of the external field [14]. Observe 
that A^z; x; C; T x ), eq. (2.5), is a special case of A^z; x; C; T x ), eq. (2.8); when T x 
reduces to T x , since, in this case V(T X ) goes to zero along with the second term in 
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(2.8). As one would expect, the external field intensity tensor is not changed by the 
transformation which takes (2.5) into (2.8), namely A^A^T^] = A^A^Tx)]. We 
emphasize that in the presence of the counterterm (2.7), as we will explicitly see, the 
correlation functions of p are completely independent of T x or T x , only depending 
on x and C. A local limit for p can then be obtained by using p, the radius of the 
circumference part of C and 5, the angular width of the region along the cut (See Fig. 
1), as regulators and at the end of all calculations, taking the limit in which p and S 
vanish. As we will see, the p correlation functions will only depend on x in this limit. 
This is the procedure we are going to adopt in order to obtain local vortex (disorder) 
correlation functions and operators. 

2.1.2) The Order Operator 

Let us introduce here the order operator which is going to be dual to p. The 
natural way of defining the order operator, which is suggested by the study of a wide 
class of models in 1 + ID [12] and also in 2 + ID [18,24], is to couple the external field 
A^ to the topological current J M in the following way 

E(x; C) = exp{ia J d 3 zJ >1 A fl } = exp{ia J d 3 ze^ A^Wp} (2.10a) 

or 

S'(x; C) = exp{ia J d 3 zW^J fl [A II }} = exp{ia J d 3 ze^W^Ap} (2.106) 

which differs from E by a boundary term. In (2.10), a is a free real parameter with 
dimension of (mass) 2. 

It turns out that we cannot obtain a local operator out of E or £' even in the limit 
where p,S — > 0. No surface renormalization counterterm will render the E correlation 
functions surface invariant in theories containing a Maxwell term in the action. This 
is so, because in contraposition to the case of p, no renormalization counterterm 
will ever be able to make the integrand in the functional integral describing E or E' 
correlation functions to depend on the external field A^ in the combination + A M , 
such that the procedure employed in [14] to show the surface invariance of p could 
be applied. The only case in which the E or E' operators could be made local would 
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be a theory containing just the Chern-Simons term in its kinetic lagrangian. In this 
case the renormalization counterterm would be a Chern-Simons term involving the 
external field. In this case however, the \x operator would no longer be local! Actually, 
as we will see, in theories containing both the Maxwell and Chern-Simons terms, we 
will have to redefine \i by adding a E' piece, in order to obtain a local disorder (vortex) 
operator. 

In order to construct a local operator oix) dual to /j,(x; C), we are going to take 
advantage of the Cauchy-Riemann equations relating the real and imaginary parts of 
the analytic function In z(see Appendix A). Let us define 

a(x°,x) = lim exp{ia [ <f£[e ij df ] arg{£ - x) + of } In |£ - x \]Wj(x°, |)]} (2.11) 

In the above expression a is a free real parameter with dimension of (mass) 5 and 
the integral is performed over the region R X (C), depicted in Fig. 3. R X (C) is the part 
of the R 2 plane external to the curve C. This consists of the arc of a circumference 
of radius p centered on x and crossing the cut of the function arg (£ — x) and two 
straight lines, along this cut, separated an angle 25 apart. 

Using the Cauchy-Riemann equation, eq.(A.2) one immediately realizes that the 

— * 

exponent in (2.11) is different from zero only on the singularities of the arg(£ — x) 
and In |<f — x\ functions, that is, the point x = (x°,x) and the cut of arg(£* — x). It is 
therefore clear that when we take out the regulators: p, 5 — > 0, a will only depend on 
x. This will be made at the end of all calculations, as in the case of p. 

As we shall see, a as given by (2.11) is the correct local operator dual to //. 

2.2) Commutation Rules. Anyon Operators 

Let us study here the various commutation rules involving the a and \x operators. 
These will enable us to construct an anyon operator out of a and /i. 
The basic equal-time commutators of the Maxwell theory are 

[W\ W j ] = [E\ E j ] = [IT, IT] = i,j = 1,2 

[W\x),Ei{y)\ = -[W\x),W{y)\ = -iS* S 2 (x - y)} (2.12) 
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Here E l = W l ° is the electric field and II 1 = — E l is the momentum canonically 
conjugate to W\ 

We show in Appendix B that the \x operator, eq. (2.6) can be cast in the form 

H&C) = exp{-ife / d 2 ^(z°, £)d ia rg(£- x)} (2.13) 

JT x (c) 

Using (2.12) we immediately see that [//,//] = [era] = 0. The \i — a commutation 
rule can be obtained by writing fi p s(x; C) = e A ( x ' c \ cr p s(y; C) = e B( - y ' c ^ and using the 
formula e A e B = e B e A e iA,B \ which is valid when [A, B] is a c-number. Making use of 
(2.13), (2.11) and (2.12), we find 

[A(x;C),B(y;C)]=ab f d\ f d 2 V d t arg^ - x)[e kl d k arg(ff - y) 

JT X (C) JRy{C) 

+di In \ff- y\]{-i)5 il 5 2 (i - ff) (2.14a) 



or 



[A(x; C), B(y; C)] = -iab / d 2 r)[e kl d k arg(fj- y) + di In \ff - y\}diarg(ff- x) 

J R y (C)nT x (C) 

(2.146) 

In Appendix C it is shown that after taking the limit p, 5 — > in R y (C), we obtain 
the following result for (2.14): (B(y) = lim Pi5 _> B(y; C)) 

IA( X - O, B(y)] = { " r ^- f > H T T f c ] (2.15) 

We therefore conclude that the equal-time commutation rule between a and ji is 
(a(y) = lim Pi<5 ^ a pS (y; C)) 

/ , p5 (x;C)a(y) = { (T ^J a; '^ eXp{ - ?27ra6 ar 9(v-*» tll'W (2.16) 

This is precisely the U(l) dual algebra appropriate for a vortex creation operator 
introduced in [14]. Taking the local limit for /i p5 , we obtain the local dual algebra 
(equal times) (fi(x) = \im p ^ n pS (x;C)) 



/i(x)a(y) = o-(y)fj J (x)exp{—i2nab arg(y — x)} (2-17) 
9 



Let us investigate now the commutation rules of a and p with the charge and 
magnetic flux operators, given respectively by 

Q = J d 2 xp(x) and $ = J d 2 xB(x) 

where 

p(x) = diE\x) and B{x) = -e^diW^x) (2.18) 

Using (2.12) we readily conclude that [a, B] = [p, p] = 0, indicating that a does 
not bear magnetic flux and p does not bear charge. The commutation rules [a, p] and 
[p,B] may be obtained by the use of the formula [e a ,/3] = e a [a,j3], valid when [a, (5] 
is a c-number. 

For e a ^ = a(x) and f3(y) = p(y), using (2.11), (2.18) and (2.12), we immediately 

get 

[a{x),/3(y)] = lim(-a) / d 2 £ [e i W i arg(i-x) + d j lim Wd^ 5 2 {£-y) (2.19a) 

Using the result (C.l) of Appendix C, and the fact that dW8{£ - y) = -&®8(£ - y), 
we obtain 

[a(x),P(y)] = 2vra 5 2 (x - y) (2.196) 



We therefore have 



or 



[°{ x )ip{y)\ = 2tt aa(x)5 2 (x - y) 
[[a(x),Q] = 2iraa(x)Q 



and 



[a(x),$]=Q (2.20) 

This result indicates that a(x) creates states bearing 27ra units of electric charge. 

Let us choose now e a ^ = p(x) and (5{y) = B(y). Using (2.13), (2.18) and (2.12), 
we get 

[a(x),P(y)\ = lim (-6) / d 2 £ Pdjargtf- £)c^ 2 (£ - y) (2.21a) 

p,6->0 JT X (C) 
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Using the result (C.2) of Appendix C, we obtain 

[a(x),(3(y)] = 2nb5 2 (x-y) (2.216) 
We therefore have, in the case of \i : 

[t,(x),B(y)] = 2nbt,(x) 5 2 (x-y) 

or 

[/jl(x), $] = 2irb/j,(x) 

and 

[ f i(x),Q] = (2.22) 

This result indicates that /i(x) creates, states bearing 2nb units of magnetic flux 
(or (-) topological charge). 

Since a bears charge and \i bears magnetic flux one should expect that the com- 
posite operator 

v( x ) = l™(r(x)n(y)z(x - y) ( 2 - 23 ) 

jj > y 

—where Z(x — y) is some renormalization factor used to absorb the short distance 
singularities in the above operator product— is an operator with generalized statistics 
or an anyon operator. Indeed, using (2.23) and (2.17) one immediately finds the 
equal-times commutation rule for ip: 

(p(x)(p(y) = (p(y)(p(x) ex.p{i2nab[arg(x-y) -arg(y-x)]} 

or 

ip(x)ip(y) = ip(y)<f(x) exp{i27iabn e[arg(x — y) — n]} (2.24) 

where e(x) = sign(x). Eq. (2.24) is the appropriate equal time commutation rule 
for anyon fields which generalizes the one-dimensional one [8,4]. Eq. (2.24) indicates 
that the spin of the anyon field ip is S v — nab, that is, S v = in agreement with 
the observations made in [3,4]. 
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Let us remark at this point that the charge bearing operator cr(x), dual to fi(x) is, 
of course, gauge noninvariant. fi(x) is essentially the generator of a gauge transfor- 
mation with parameter arg(y — x), as can be clearly seen from (2.17) and (2.16). A 
gauge invariant operator would commute with n(x) and therefore could never possibly 
be its dual! Indeed, starting from (2.11) and using (C.l) we see that under the gauge 
transformation — > + A, we have a(x) — > exp{i27raA(x)}(x(:r). With respect 
to gauge transformations and charge bearing, the order operator cr(x) is very much 
like the scalar field <fi of the abelian Higgs model or the electron field of Quantum 
Electrodynamics. 

2.3) Disorder (Vortex) Correlation Function 

Let us evaluate here the two-point correlation function of fi within the euclidean 
functional integral framework. Going to the euclidean space (ix° — > x 3 ,iW° — > 
W 3 ,iA° -> A 3 ), we will have 

iS = ij d 3 z{- l -W^W^} - -J d 3 z E ^W, v W^ = 

- J d 3 z E ^W^W u = -S E (2.25a) 
where P^ v = —U E 5^ V + & l d v ' . The exponent of fx in(2.6) will transform as 

~ J d 3 zW^A, u - / d 3 z E W^A liV = -bj d 3 z E W,P^A v (2.256) 

We will also use a gauge fixing term of the form 

iS GF = i J d 3 z[-^(d,Wn 2 } - ~\ J d 3 z E (d,W») 2 = 

~ J d 3 z E W^W v = -S e ,gf (2.25c) 

where G^ u = —£d^d u and £ is an arbitrary real parameter. The renormalization 
counterterm (2.7) transforms as 

iS R = iJ d 3 z[-^A^A^] - -J d 3 z E jA liU A^ = 

- J d 3 zJ^A,P^A u = -S e ,r (2.25d) 
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Using (2.25) we can write the following expression for the euclidean two-point 
correlation function < ji{x) ji* (y) > (we henceforth will drop the subscript E): 



\w,G^W u + W, u A^(z; x, y) + ±A%,(z; x, y)}} (2.26a) 

< fi(x)fj*(y) >= lim Z- 1 [ ZW„exp{- f d 3 '' ziw^P^ + G» V ]W U 
p,s^o J J 2 

+W,P^A u (z; x, y) + \a\ v {z- x, y)}} (2.266) 

In these expressions, Z is the vacuum functional and A^z; x, y) = b[A^(z; x; C; T x ) — 
A^z; y; C; T y )] (the minus sign corresponds to the fact that we have the adjoint op- 
erator n*(y)). The mixed — A^ terms in the exponent in (2.26) come from the 
fj,(x) fj,* (y) operators and the last term is the surface renormalization counterterm. 

Observe that, as we put forward above, the integrand in (2.26) depends on the 
external field in the combination + A^. Under the change of variable — > 
W 7 ^ + dpid, with oo given by (2.9), we can arbitrarily change the surfaces in A^ [14]. 
The gauge fixing term would be invariant under this transformation, provided we add 
the ghost term to (2.26) [21]. For simplicity reasons, however, we will neglect the 
ghost term in the abelian theories considered here. The functional integral in (2.26) 
can be evaluated with the help of the euclidean correlation function < W^W V >= 
[P^ u + G^Y 1 . This is given, in momentum space, by 

< W»(k)W v (-k) >= I ' + ; P^(k) = 5^k 2 - k»k u (2.27a) 

k £ k 

and in coordinate space by 

< WHx)W v (y) >= lim [P** - -<W][— - ^ ~ V \ ; P^ = -5^ + d»d v (2.276) 
m->o £ m 8tv 

Here m is an infrared regulator used to control the singularities of the inverse 
Fourier transform of (see(4.15)). 

Using the results of Appendix B, we conclude that each piece of the second term 
in the exponent in (2.26b) can be put in the form 

-b J d 3 zW^ v A u = ~J dhW^d^ = J d 3 zB^W\ (2.28a) 
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where 

B,(z; x) = bf d 2 i a df arg{i- x)F a ? ,5 3 (z-Z), (2.28b) 

in which 

F afB 7 = d a 5 p 7 - d p 5 a 7 
Performing the functional integral in (2.26), we get 

< fi(x)fi*(y) >= lim exp{J f d 3 zd 3 z'B^(z;x,y)B u ( y z';x,y) 

x [P^ + G^]-\z ~ A ~ Sr[A^}} (2.29) 

where B^z^x^y) = B^(z]x) — B^z^y) and Si?L4^] is the last term in the exponent 
in (2.26). Using (2.28b), (2.27b) and integrating over z and z', we obtain 

<n(x)(i*(y) >= Jim exp{^- J2 e * e i / rf2 ^ / d 2 r] a dl arg(^-Xi)d { ^ arg(rf- 
M^-o 2 iJ=l Jt x . (c) Jt Xj(C) 

* F & ^5) ^ P (0 ~ £ - ^ - ^1) (2- 3 °) 

In the above expression, x± = x, e± = 1 and rr 2 = y, e 2 = — 1. 

We immediately see that, since F a/3 x d x = 0, the gauge dependent part of < 
W^W U > gives no contribution to < pL/i* >, as one should expect for a gauge invariant 
operator as //. The only contribution comes from the gauge independent part of 
< W^Wy >, which is proportional to —U5^ v . 

Let us now make use of the very useful identity 

- DF (0 ° F (v) ^ = -□^V^ ) %^ - 4 } %) (2.31) 

which can be readily demonstrated by expanding the e' s in terms of Kronecker 
5's. Inserting (2.31) in (2.30) and observing that 

_□[! - M] = -OF-^l] = jr-i[J_] = i im 1 (2.32) 

(here e is an ultraviolet regulator used to control the short distance singularities of 
the inverse Fourier transform of p-) we get 

b 2 2 

< fi(x)fi*(y) >= lim exp{-— 
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/ d% I d 2 Va d^arg(£- x^ argtf - Xj) 

jT Xi(C) jT xj{C) 



Here we used the fact that d^F(^ — 77) = —d(£\F(£ — rj). 

Observing that lim e ^ ~~ □ [47r [ | rc | 2 + | e | 2 ]] ~ i = 5 3 (x) and using the results of Ap- 
pendix B, we see that the first term in the exponent in (2.33) can be written as 

% E / d% I d 2 Va d^arg(i- x^dfargtf- £,)e^V^ 3 (£ - r,) 

i j = i-JT Xi (C) JT X .{C) 

= - A j dhdh'^A^x^y) - (pt <- z/)]5 3 (^ - ^[dMV^y) - (/x <- !/)] 

= J / <FzA%,(z; x, y) = S R [A,} (2.34) 

We immediately conclude that this term is canceled by the surface renormalization 
counterterm in (2.33). 

The second term in the exponent in (2.33) can be evaluated by making use of the 
results of Appendix C. It is surface independent and, in the limit p — > only depends 
on the points x and y. Indeed, according to (C.2) and already using (2.34) 

b 2 2 n 
< fi(x)n*(y) >= limexp{-- ]T ei6j .[ — ———-— ]} (2.35a) 

i,j=l I 1 3 1 ' I I J 

or 

< n(x)n*(y) >= limexp{7r& 2 [ 1 - -]} (2.356) 

e * 1*^ y 1 ^ 

We now see explicitly that the renormalization term (2.7) indeed, renders the pi- 
correlation function surface independent. 

The ultraviolet singularity in (2.35) may be eliminated by a renormalization of the 
disorder (vortex) operator: 

Tib 2 

Hr{x) = hm p(x) ex P[^] ( 2 - 36 ) 
We therefore finally arrive at the result 

7Tb 2 

< VrWvUv) >= exp{^— } (2.37) 

\x y\ 
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for the disorder (vortex) euclidean correlation function in the Maxwell theory. Observe 
that 

lim </i fl (iK(»)>=l (2.38) 

\x-y\-KX 

This indicates that hr does not create genuine vortex excitations, in the true Maxwell 
theory. One should expect this result since already at the classical level the theory 
does not possess topological excitations. Another way of describing this fact is to say 
that there is a vortex condensate in the Maxwell theory. In section 3 we will examine 
how this can be changed in the Maxwell-Higgs theory. 

2.4) Order Correlation Function 

Let us evaluate now the two-point euclidean correlation function of a. From(2.11) 
we see that there is no change in a as we go to euclidean space. Using the results of 
the previous subsection for the euclidean action S and gauge fixing term Sqf, as well 
as the expression of a(x) eq.(2.11), we may write 

<a(x)a*(y) >= lim Z' 1 f DW^exp{- f d 3 z^-W^P*" + G^]W V + C^z; x, y)W»}} 

(2.39) 

In this expression, C^(z;x,y) = [C^z; x; C; T x ) — C^(z;y;C;T y )](as in the case 
of //, the minus sign corresponds to the fact that we have the adjoint operatorcx*(|/)) 
where 

C^{z] x] C; T x ) = 

= { ™f Rx (C) d 2 i[e k ^arg(i- x) + dP In |£ - x \}5 3 (z - £) for p = i 

\ for fj, = 1 ' ' 

Performing the functional integral in (2.39), we get 

< a(x)a*(y) >= lim exp{i / d 3 zd 3 z'C,(z- x, y)C v {z'- x, y)[P^ + G^]-\z - z')} 

(2.41) 

The a correlation function, of course is not gauge invariant since a itself is not. 
In order to obtain the gauge independent part of < a a* >, let us use in (2.41) just 
the gauge independent part of < W^W^ >, namely, the 5^ v - proportional part of 
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< W^W P > in (2.27b), which is explicitly given by (2.32) (times 5^). Inserting this 
and (2.40) in (2.41), and integrating on z and z' , we get 



< a(x)a*(y) >= lim exp{- — V e^e,- / d 2 £ [ d 2 nx 
KJ yy> M~o i X 2 1 3 jR Xi{ c) ^Jr X] (c) 1 



[cPdfpargg-fy+dfa In [e^d^arg(ff-x J )+cP {v) In |q-s 3 -|] [ _ ^ + ]} 

(2.42) 



Let us now use the identity 



e -m\\ 



V ri i2. i 1211/2 -^1 ]+ 



[-_J__ + m ] e-»0-l-^ (2.43a) 



M 2 + |e| 2 ] L [|z| 2 + |e| 2 ]V 
which in the limit m — > reduces to 

9 ^[__L + [|af + | c |2]l/2] + * * (2.436) 



[|x| 2 + |e| 2 ]V2 ttt. J J [|x| 2 + |e| 2 ] 3 / 2 

Insertion of (2.43b) in (2.42) produces two pieces. The one coming from the second 
term in (2.43b) vanishes in the limit p, 8 — > as is shown in Appendix C. The one 
coming from the first term in (2.43b) can be computed with the help of (C.l), after 
using the fact that dfoF(£ — rj) — — c^F(£ — 77). We therefore conclude, according 
to (C.l) that 

< a(x)a*(y) >= lim exp{^ £ e ^( 2 7r) 2 [-i + [\ Xi - x 3 \ 2 + |e| 2 ] 1/2 ]} (2.44a) 

m,e^0 S7T • 777. 

or 

< cr(x)cr*(y) >= exp{-7ra 2 |:r - y|} (2.446) 

Observe that the infrared singularities at m — > completely cancel. If we were 
calculating a charge nonconserving correlation function (as < era >, for instance) we 
would have = and the — factors would no longer cancel, implying < aa >= 
as m — > 0. We see that, as in 1 + ID [4], an infrared singularity is responsible for 
enforcing the selection rule for the dual operators. Observe that 

lim < <r(x)cr*(y) >= (2.45) 

|rr-2/|->0 
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indicating that the charged states created by a are orthogonal to the vacuum and 
charge is a conserved quantity as one would expect in the pure Maxwell Theory. 

2.5) The Mixed and Anyon Correlation Functions 

Let us examine now the mixed order- disorder euclidean correlation function from 
which we will be able to obtain the anyon correlation function. 
Combining (2.26) with (2.39), we can write 

< a(x 1 ) f i(x 2 )a*(y 1 ) f i*(y 2 ) >= lim Z' 1 f DW^exp{- f d 3 z\]-W,[P^ + G^]W V + 

p,o^0 J J Z 

C^z; x u Vl ) W* + W^A v {z- x 2 , y 2 ) + \a\ v {z; x 2 , y 2 )}} (2.46) 
The functional integral can be evaluated as before, yielding 

< a(x 1 )fJL(x 2 )<T(y 1 )fi(y 2 ) >= lim exp{J f d 3 zd 3 z'[B^(z; x 2 ,y 2 ) + C^z; x 1 ,y 1 )] 

p,o^0 Z J 

x [B u (z' ; x 2 , y 2 ) + C„(z'; x 1 ,y 1 )} [P^ + G^]~\z - z') - S R [A,(x 2 , y 2 )}} (2.47) 

The BPB and CPC terms were evaluated before. The novel terms which appear 
here are the BPC terms(already considering the two of them): 

r 2 
BPC = / d 3 zd 3 z'B^{z;x 2 ,y 2 )C u (z';x 1 ,y 1 )[P^ + G» u ]-\z - z') = lim iab V 

J ' p,S,m^0 . j=1 

x / d 2 ^ I d 2 rid^ ] arg(i- f^d^ arg{ff - s,) + d i In \rj - sA] 

JTrAC) JRsAC) 



xF (0 ^ P ii) - -c d h d yi- - ^7^] k k = 2 ) ( 2 - 48 ) 



^C) JR S] {C) 

Here r\ = x 2 , r 2 = y 2 , Si = X\ and s 2 = y\. As before e\ = 1 and e 2 = —1. As in 
the case of < jiji* >, the gauge dependent part of the propagator in (2.48) gives no 
contribution, because F^ u a d a = 0. Only the — 05 m part contributes: 

F pu ^ P ai = _ U (qp S * _ grs^J 7 - 1 ^] = {jpsr - d v 5^)T~ l [^\ (2.49) 

/C K 

The t^-part of (2.49) vanishes because d 2 ^ is orthogonal to R Sj (C)- The remaining 
part gives 



2 . 

BPC = lim iab Y e t e, / [e^d^arg 



(j7-Sj) + 9(„) ln|?7-s 
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^kic/^^-^^W^ (2 - 50) 
The above integrals are evaluated in Appendix D. Using (D.l), we get 

2 

BPC = iirab €i€jarg(sj — fj) 

or 

-BPC = inablargtx! - x 2 ) + arg(y x - y 2 ) - arg{x x - y 2 ) - arg(y x - x 2 )\ (2.51) 



Combining (2.51) with the previous results for the BPB and CPC terms, we obtain 

irb 2 

< <T{x 1 )fj, R (x2)<T*(y 1 )fi* R (y 2 ) >= exp{-7ra 2 |a;i - y ± \ + 

\x 2 -y 2 \ 

+inab[arg(x 1 - x 2 ) + arg(y x - y 2 ) - arg(x x - y 2 ) - arg(y x - x 2 )]} (2.52) 



The anyon correlation function can now be obtained in a straightforward manner. 
Introducing the anyon field 



ip(x) = lim_ ai^XijiiR^x^ exp[— inab arg{x\ — x 2 )\ (2.53) 

X\ — >X'2 = X 

we immediately get 

nb 2 

< Lp(x)Lp*(y) >= exp{— na 2 \x — y\ + — inab[arg(x — y) + arg(y — x)]} (2.54) 

\x y\ 

Observe that < (pip* > is multivalued, the ambiguous phase being exp[i27r(7ra6)]. As 
in the case of 1 + IP [4], we interpret this fact as being the manifestation in the 
framework of euclidean correlation functions of the nontrivial commutation rule of 
ip. Of course, the same functional integral describes the two possible orderings of 
operators in the left hand side of (2.54). Making the analytic continuation from 
euclidean to Minkowski space from each sheet of (2.54) would reproduce each possible 
ordering of </?[4]. This, of course implies that the spin of ip is s v = nab, confirming 
the result we found above by direct computation of the ip commutator. 

The (/^-correlation function (2.54) decays, at large distances as exp[— ira 2 \x — y\\. 
This implies the anyon states created by ip possess a mass M = na 2 . 
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3) The Maxwell-Higgs Theory 
(Abelian Higgs Model) 

3.1) Introduction 

Let us study in this section the properties of the dual operators a and \i introduced 
before, in the framework of the Abelian Higgs Model (AHM), described by 

C = -\wl„ + |ZVI 2 - ffVV - \^<P? (3.1) 

where = + ieW^. 

The theory exists in two phases: i) a symmetric phase, for g 2 > 0, where < <fi >= 
and is massless and ii) a broken phase, for g 2 < 0, where < <fi >= k ^ and 
acquires a mass M = en through the Higgs mechanism. 

Let us introduce in this theory the a and \x operators defined in section 2. Since the 
equal time commutators of the AHM are exactly the same as in (2.12), we conclude 
that all commutation rules (equal times) evaluated in section 2.2 for the pure Maxwell 
theory remain valid here. 

As the surface independence of \x is concerned, let us note that the same countert- 
erm (2.7) will be sufficient to make the /^-correlation functions surface independent 
here, because also in the AHM it will make the integrand in the functional integral 
defining < fifi* > to depend on the external field through the combination + A^. 
The remaining three last terms in (3.1) will be invariant under the change of variable 
(involving cu, eq. (2.9)) needed to show surface independence. As was shown in [14], 
/i is the operator which creates the quantum states associated to the classical soliton 
(vortex) solution of Nielsen and Olesen [22]. 

In this interacting theory, of course, we will no longer be able to compute exact 
correlation functions. Instead, we will evaluate the long distance behavior of them. 
As we will see, this will be enough to obtain interesting physical consequences. 

3.2) The Unbroken Phase 

Let us study here the symmetric phase (g 2 > 0), in which < >= 0. We start 
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with the fx correlation function which is given by an expression similar to (2.26): 
<fx(x)fi*(y) >= lim Z- 1 f DW^exp{- f dh^W^ + G^]W V + \D^\ 2 + V(<f>) + 

W^A v (z; x, y) + \a\ v {z- x, y)}} (3.2) 

Using (2.28), we see that < fx/x* > is obtained by coupling the external field 
B^(z; x, y) to the AHM lagrangian, in the way given by (2.28a). We therefore conclude 
that 

< fx(x)fx*(y) >= exp{A(x, y) - S R [A^x, y)}} (3.3) 

where, using a diagrammatic language, A(x,y) is given by the sum of all Feynman 
diagrams with the field B^{z\x,y) in the external legs. 

We are only interested here in the long distance behavior of < fx/x* >. As was 
shown in [21], only two-leg graphs contribute to (3.3) in this limit. We are going 
to evaluate A(x,y) by making a loop expansion (in powers of H). The only vertex 
involving the external field is given in Fig. 4a. The lowest order (two-leg) graph 
contributing to A(x, y) is given in Fig. 4b. Inserting the euclidean propagator for W^, 
eq. (2.27) we immediately see that the graph of Fig. 4b. is identical to the first term 
in the exponent in (2.29). It follows that in this order of approximation (0(|), the 
large distance behavior of < fifi* > is given exactly by the same expression as in the 
Maxwell theory, namely 

< /x*(*K(y) > l ^°° exp{^} **?tr 1 (3.4) 

\x y\ 

We have renormalized fx in the same way as in (2.36). 

From (3.4) we conclude that in the symmetric phase of the AHM, in analogous 
way as in the true Maxwell theory, the disorder operator fx does not create states 
orthogonal to the vacuum, that is, genuine excitations, as one should expect. Also 
the symmetric phase of the AHM can be viewed as a vortex condensate. 

Let us turn now to the order correlation function. This is now given by 

< a(x)a*(y) >= lim Z' 1 f DW^exp{- f d 3 [ z[-W M [P^ + G^]W U + 
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+ |ZV>| 2 + V(<f>) + C,(z; x, y)W»}} (3.5) 
where C^z^x^y) is defined in (2.39) and (2.40). Now, 

<a(x)a*(y)>=exp{T(x,y)} (3.6) 

where T(x, y) is the sum of all Feynman graphs containing the field C^(z; x, y) in the 
external legs. Again, only two-leg graphs contribute to the long distance behavior of 

< (jo-* > [21]. The only vertex involving C M is the same as for and is shown in 
Fig. 4a. Again, making an expansion in loops, the lowest order graph will be the one 
of Fig. 4b. Inserting the W^-propagator, eq.(2.27), we conclude that in this order 
(O(jr)) the large distance behavior of the gauge invariant part of < era* > is given by 
the same expression as in the Maxwell theory, namely 

. / \ */ \ ^ \x-y\-HX> r 2| n \x-y\->oo , , 

< a(x)a (y) > — > exp{— na \x — y\\ — > (3.7) 

The long distance behavior of the anyon field could be obtained as well, by ex- 
changing B^iz-.x.y) for B^{z]x,y) + C^z\x,y) in (3.2). In the lowest order of ap- 
proximation, we would get < ip(x)ip*(y) > behaving asymptotically as (2.54). 

3.3) The Broken Phase 

3.3.1) The Order Correlation Function 

Let us consider now the case in which g 2 < and < >= k = (^M-)^. We 
are going to write in terms of its real components: = ^{<f>i + ifo) and choose 

< 0i >= k and < 02 >= 0. 

Let us compute here the order correlation function< aa* >. This is given by 
(3.5), and (3.6). Again, in lowest order < aa* > will be given by the graph of Fig. 
4a. There is, however, an important difference. We must shift the Higgs field around 
the vacuum value. This will generate the following quadractic part for the euclidean 
lagrangian C[W^ 0i, <j> 2 ] 

£ (2) [^,0i,0 2 ] = + l -M 2 W,W, + i(^0x) 2 + \m\^\ + \{d^ 2 f + MW^th 

(3.8) 
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where M = en and mi = 2\g 2 \. 

Using the Lorentz gauge (£ — > oo), we see that PU M and 2 decouple in the quadrac- 
tic part (since d^W^ = 0). The euclidean propagators for W^, 4>\ and 2 are now 
given respectively by 

D "- {k) = & [ Fp + W) + ¥WTW) ] (3 ' 9tt) 

Inserting the W^-propagator in the graph of Fig. 4b, we see that, the lowest order 
(0(|)) contribution to large distance behavior of the gauge invariant part of < oo* > 
is given by an expression identical to (2.41) and (2.42), except for the fact that the last 
term between square brackets in (2.42), namely is exchanged by ^ r ~ 1 [ fc2+ 1 M 2 ], 

that is 

1 e -N 2 +N 2 ] 1/2 
^ v r- l [- — ] = lim 5^ — --^ (3.10) 

l k 2 + M 2> e^O [| x |2 + | e |2]l/2 V J 

Using the identity (2.43a) and following the same steps which led us to (2.44), we 
immediately conclude that in lowest order (0(^)), 

< o(x)o*(y) > l ^°° exp{^[e- M ^l - 1]} (3.11a) 

or 

< o R (x)o* R (y) > l ^°° exp{^ e - M l^^} 1 (3.116) 

where o R = aexp[ff^]. 

Observe that in the broken phase, the charge bearing operator o no longer cre- 
ates states orthogonal to the vacuum. This corresponds to the screening of charge 
associated with the mass generation to the field through the Higgs mechanism. 
Notice that in the limit M — > 0, < oo* > reduces to the expression found in (3.7) for 
the symmetric phase and the zero mass singularities cancel in charge selection rule 
respecting correlation function. 

3.3.2)The Disorder (Vortex) Correlation Function 
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Let us evaluate here the /i-correlation function in the broken phase of the AHM. 
As we infer from (3.2), we can write < /i/i* > as 

< n(x)»*(y) >= lim Z- 1 [DW.expi-SlW^ + A^^D^-SGFlW,]} (3.12) 

where SpH 7 ^, 0, D^ifi] is the action associated with (3.1). 

Performing the change of variable — > + A^(z\ x, y), we get 

< n(x)n*(y) >= lim / D^expi-SiW^^D^] - S GF [W, - A,}} (3.13) 

where = d^+ie[W^ — A^(x, y)]. It will be more convenient to use (3.13) to compute 
< /i/i* > in the broken phase. Now, we can write 

< n(x)n*(y) >= exp{A(x,y)} (3.14) 

where A(x, y) is the sum of all Feynman graphs with the A^(z; x, y) field in the external 
legs and computed with the Feynman rules coming from (3.13). We are again going 
to make an expansion in loops and furthermore, in powers of M, the gauge field mass. 
Before evaluating A(x,y), we shift the Higgs field around the vacuum value in (3.13). 
This operation commutes with the change of variable — > + A^. The vertices 
relevant for the lowest order computation of A(x,y) are depicted in Fig. 5. Their 
contribution, in this order, is given by the graphs of Fig. 6. The gauge dependent 
terms, which appear in Fig. 6a, vanish for all values of £ as can be easily seen from 
(3.9a). The first nonzero contribution for A(x,y), therefore is given by the graphs of 
Fig. 6b, which are of the order O(^). Using (3.9), in the graphs of Fig. 6b, we find 

A(x,y) ~ J d 3 zdh'A,(z ] x,y)A u (z' ] x,y)P^F(z - z>) (3.15) 

where F(z - z') = F' 1 ^} = A 2 (z - z') and is given by (2.32). Eq. (3.15) can be 
written as 

A(x,y) l ^ ° ^ J Szd i z'e^ x d^A^z ] x,y)F{z-z')e^ x df ) A u {z' ] x ) y) 

Using the result of Appendix B, we get, after integrating over z and z' and using 
(2.32): 

~ \ x - y \^oo M 2 b 2 2 



^|zj/|_^oo I' [ d 2 ri a d^parq(f—Xi)d^arq(ff 

M^o 2 ^ 1 3 Jt x .(C) ^JT Xj (C) la v ^ 1,(5 yW 
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x^^lir^ 12 i i 211/2 ] (3-16) 

47r[|£ - rj\ 2 + |e| 2 ]V2 J v 7 

Making use of the identity (2.43b) we see that (3.16) contains two terms, corre- 
sponding to the two pieces in the right hand side of (2.43b). The second one is surface 
dependent but, as was shown in [17], for i ^ j, it vanishes in the limit \x — y\ — > oo. 
This therefore does not contribute to (3.16). The i — j term is just a self-energy that 
will renormalize \i [17]. The first term of (3.16) can be evaluated immediately with 
the help of result (C.2) of Appendix C: 

A'^-lim^ ± ^[-i + ^-^f + Hr 2 ] (3.17a) 



or 



A(x,y) lx ^°° -nM 2 b 2 \x-y\ (3.176) 



We therefore conclude that 

< ii R {x)ii%{y) > l ^°° exp{-7rM 2 6 2 |x - y\} (3.18) 

We see that in the broken phase the vortex operator fx indeed creates states orthog- 
onal to the vacuum, i.e., genuine excitations in the spectrum. According to (3.18), 
the mass of these excitations is M v = irM 2 b 2 , up to the order of approximation we 
are working (0(^=p)). It is interesting to see that choosing b 2 = e~ 2 [21], the mass of 
the quantum vortices created by /j, coincides with the classical vortex energy found in 
[23]. In ref. [21] one-loop corrections to the above result were also evaluated. 

Our study of the AHM and Maxwell theory clearly exposes the reason why we call a 
and \i "order" and "disorder" operators, a has a nonzero vacuum expectation value in 
the broken (ordered) phase, while the vacuum expectation value of \i is different from 
zero in the symmetric (disordered) phases and vice-versa. In the broken (ordered) 
phase fj, creates genuine quantum soliton states. The transition to the symmetric 
(disordered) phases can be viewed as a condensation of topological charge in the same 
way as the transition to a broken (symmetric) phase can be viewed as a condensation 
of charge as occurs in superconductivity. We see that in a certain sense, the AHM 
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is quite similar to the Ising model, the phase transition which generates mass to the 
gauge field being analogous to the phase transition which takes place in the latter. 

A final word is in order about the procedure adopted for the computation of 
< /ipL* > in the broken phase of the AHM. Observe that we shifted around 
before evaluating this correlation function, contrary to what we did in other cases. We 
proceeded this way because without shifting, the lowest order in our approximation 
scheme would produce a trivial result as the computations would be effectively done 
in the Proca theory. We would have to go to higher orders in order to get a nontrivial 
result, while with the shifting trick we already get it in lowest order. Of course the 
two ways of calculating should produce the same result if one was able to sum the 
whole perturbation series. 

4) The Maxwell-Chern-Simons Theory 

4.1) Order and Disorder (Vortex) Operators 

Let us study now the Maxwell-Chern-Simons (MCS) theory, described by 

C = —W^W^ + 6 -e^ a W,d a W p (4.1) 
which leads to the field equation 

d v W^ v = 9e^ a(3 d a W,_ 

The parameter 9 has dimension of mass. 

In addition to the identically conserved current (2.3), let us introduce the two 
indexes current Jq V which generalizes (2.4): 

Jf = W» v - 9e> M/a W a (4.3) 

Jq U is conserved as a consequence of the field equation (4.2) in analogy with (2.4). 

Let us introduce also the disorder operator pig which generalizes (2.6) for 9 ^ 0, 
by coupling Jq V to A^: 

ib f ~ 
Ho(x; C) =exp{-— J d 3 zJg l ' A^} 
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(4.2) 



or 

He(x; C) = exp{| f SzW^A^ + ib9 J dhe^W^Ap} (4.4) 
We see that fig = fi(b)Tl (9b) , where E' is given by (2.10b). 

Making use of the surface renormalization counterterm (2.7), we immediately see 
that the integrand in the functional integral defining correlation functions of fig will 
depend on through the combination + A^ and therefore, we expect these cor- 
relation functions to be surface independent as before. We are going to see explicitly 
that this is indeed the case. 

The a operator needs no modification in the MCS theory. It is given, as before, 
by eq. (2.11). 

4.2) Commutation Rules 

Let us obtain here the relevant commutators involving the a and fig operators in 
the MCS theory. 

The momentum canonically conjugate to W % is now given by 11* = — E % + |e y 'W J ', 
with E % = W l ° [1]. The basic commutators of the MCS theory are [1] (equal times) 

[W\ W j ] = [IT, IT] = 

[W i (x),E^y)} = -[W\x),W (y)} = -iff* (2 -ft (4.5) 

[E\x),E^y)] = -t6e^5 2 (x-y) 

Observe that contrary to the Maxwell case, the electric field E % no longer commutes 
with itself. 

Using the results of Appendix B, we see that fig can be written as 

He(x; C) = exp{-ib f d 2 ^[E\x°, |) + Q^W 3 (x\ £)}d ia rg(£- x)} (4.6) 

Using (4.5) we immediately see that the commutation rule between a and fig in the 
MCS theory is identical to that in the Maxwell theory involving a and fi, namely 
(2.16) or (2.17). Also in the MCS theory, we have [a, a] = 0. fig, however, no longer 
commutes with itself. Writing fig(x; C) = e A ^ x ' ,c \ we have according to (4.5), 

[A(x; C), A(y; C)\ = lim -i6b 2 f d 2 £ [ d 2 ^ arg{i - x)df ) arg{ff- y) 
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xe ij 5 2 {£-ri) (4.7a) 

or 

[A(x; C), A(y; C)] = lim -16b 2 [ d 2 ^d t arg(£- x^arg^- y) (4.76) 

p,5^0 JT x (C)nT y {C) 

It is shown in Appendix C, that, after taking the limit p, 8 — > 0, in T X (C) and T y (C), 
eq. (4.7b) is given by 

[A(x; C), A(y; C)\ lim —i27c9b 2 [arg(x — y) — arg(y — x)] 

= -i2n(7v6b 2 ) e[arg(x - y) - n] (4.8) 
where e(x) = sign(x). It follows that, at equal times 

He{x)ne{y) = iie(y)M x ) exp{-i2n(n9b 2 )e[arg(x - y) - ir]} (4.9) 

indicating that in the MCS theory the disorder (vortex) operator is itself anyonic, 
carrying statistics = n9b 2 . 

Let us investigate now whether the commutation rules of a and fie with the charge 
and magnetic flux operators, (2.18), are changed in the MCS theory. 

According to (4.5), we immediately see that the commutation rules involving a 
are the same as in Maxwell theory, namely, eq. (2.20) remains valid, indicating that 
a bears charge but not magnetic flux. The commutator between fig and the magnetic 
flux (topological charge) operator is also identical to the one in Maxwell theory. Let us 
evaluate the commutator between fi e and the charge operator. Writing fig(x) = e a ^ 
and p(y) = djE j (y) = (3(y), we have, according to (4.5) 

[a(x),P(y)\ = lim -ib f d 2 £ df ] arg{(- x)df [-i9e ij + i9e ik 5 jk ]5 2 ^ - y) = 

(4.10) 

We see that also in MCS theory we have \jj,g, Q] = and eq. (2.22) remains fully 
valid. 

A very interesting phenomenon occurs in the MCS theory, involving the statistics 
of fig. In spite of bearing only magnetic flux (and no net charge) fig has generalized 
statistics S^ e = n9b 2 . This can be understood as follows. In the pure Maxwell 
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theory, the two quantities determining the statistics were the magnetic flux and the 
charge, whose densities were given respectively by the zeroth components of — J 71 and 
J M = d u J Ufl , where J Ufl is given by (2.4). In the MCS theory, however, must be 
exchanged by , eq. (4.3). It follows that in addition to the magnetic flux, the 
other quantity relevant for determining the statistics in the MCS theory will have its 
density given by the zeroth component of Jg = d^Jg^, or Jg = p — 9B. In other 
words, the statistics in the MCS theory will be determined by the product of $ and 
the effective charge Q = Q — Indeed, we see that according to (2.22) (which is 
also valid for fig) fig carries 2-nb units of $ and —2it9b units of Q, implying that it will 
have statistics S^ e = = ^Ob 2 (Observe also the sign in (4.9). 

A composite anyon operator can also be constructed in the MCS theory through 
an expression like (2.23) or (2.53). Using the a — fie and fig — fig commutation rules, 
we immediately find the equal-time commutator 

ip(x)ip(y) = f(y)ip(x) exp{i27i[iiab — n9b 2 ] e[arg(x — y) — it] } (4. 1 1) 

This indicates that the composite anyon field has now statistics = 7rb\a — 9b\. The 
field if carries magnetic flux $ = 2irb and effective charge Q = 2n(a — 9b). We see 
that here also the formula S v = holds true. For the special case a = 9b, ip will 
be a boson, corresponding to the fact that Q = for this value of a. 

4.3) Order Correlation Functions 

Let us compute now the euclidean order correlation function in the MCS theory. 
In order to do this we will need the analytic continuation of the Chern-Simons action 
to euclidean space: 

iScs = y/ dhef^W^Wu - - / d 3 z E ^-e^W,d a W u 

= - e -fd 3 z E W^W u = -S C s,E ; C^ = -ie^d a (4.12) 

Following the same steps as in the case of Maxwell theory and using (4.12), we arrive 
at an expression similar to (2.39) for < go* >: 

< a(x)a*(y) >= lim Z' 1 f DW^exp{- f d^z^-W^P^ + 9C^ + G*"]W V + 
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C„(z;x,y)W»]} (4.13) 

where Cn{z;x,y) was defined in (2.39-40). 
Integrating over we readily obtain 

<a(x)a*(y) >= lim exp{i f <^z<^^CJz-,x,y)C v (z'-,x,y)[P' M/ +eC' u/ +G' u/ ]- 1 (z-z')} 
p,s—*o 2 J 

(4.14) 

where [P^ + <pC» u + =< WW > M cs is the euclidean propagator of the 

field in the MCS theory. This is given in momentum space by 

< W*(k)W(-k) >mcs= [P^(k) - Oe^K] [ k2{k2 \ 02) ] + ^ (4.15a) 
and in coordinate space by 

< W(x)W(y) > MCS = [P^ + ^e^d a ][ ]'J % - lim -d^[- - ^-^] 

(4.156) 

Here, the second expression between brackets is ^ r ~ 1 [ fc2 ( fc 2 +6> 2) ] and the third is •T 7-1 ^] 
(as before, m is an infrared regulator). Observe that (4.15) reduces to (2.27) in the 
limit 6 — > 0(m = in this limit). 

As before, we want to extract the gauge independent part of < era* >. This will 
be achieved by inserting the gauge independent part of < W*W U >mcs, namely the 
5^ u and e ^ au - proportional terms of (4.15), in (4.14). The 5^- proportional part of 
(4.15) is 

1 _ p -6\x-y\ i i -6\x-y\ 

-□[- — ] = -DF- 1 \ ]=J r ~ 1 \ — - — 1 = — (4 16) 

[ Aix6 2 \x-yY [ k 2 (k 2 + 6 2 y [ k 2 + 6 2i 47r\x-y\ 1 ' ; 

We can now use the identity (2.43a) and go through the same steps which led to 
(2.44). As before, the second term in (2.43a) will give a null contribution in the limit 
p, 5 — > 0, as is shown in Appendix C. The e^^-proportional term in (4.15) also gives 
a vanishing contribution in this limit, for the same reason (see Appendix C). 

The contribution from the first term in (2.43a) can be evaluated in the same way 
as we did in (2.44), leading to 

< a(x)o*(y) >= limexp{^ £ e^Tr) 2 ^"*^ ~ X ^ + |e|2]} ]} (4.17a) 
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2 

< a{x)a*(y) >= exp^te^^ - 1]} (4.176) 

9 

or 

< a R (x)o> R (y) >= exp{^e-^^l} (4.17c) 

where <jr = crexpf 2 ^]. Observe that (4.17) reduces to (2.44) in the limit 9 — > 0. 
Now, 

lim < (J^K^) >= 1 (4.18) 

|X-J/|-KX> 

This result reflects the screening of charge associated with the mass generation to the 
gauge field induced by the Chern-Simons term. 

4.4) Disorder (Vortex) Correlation Function 

Let us evaluate here the \iq- correlation function in the MCS theory. Using (4.4), 
(4.12), (2.25) and (2.7), we can write 

< ne{x)n* e {y) >= lim Z' 1 f DW,exp{- f Sz[\w,[P^ + 9C^ + G^]W V + 

W^Pi" + QC^\A v {z- x, y) + \a\ v {z; x, y)]} (4.19) 
where A^z^x, y) was defined in (2.26). 

Notice that, since Scs[An] = because A^ only has the \x = 3 component different 
from zero, the integrand in (4.9) depends on A^ through the combination W^+A^. The 
method employed in [14], therefore, can also be used to prove the surface invariance 
of < noul > in the MCS theory. 

Integrating over W M with the help of (4.15) and using the results of Appendix B, 
we get 

< ne{x)n%(y) >= lim exp{^ j zd z z'\B^z; x,y) + D^z; x,y)} 

x [B u (z'; x, y) + D u (z'; x, y)\ [P^ + 9C^ + G^\-\z ~ A - Sr[A,]} (4.20) 

In the above expression, B^(z;x,y) was defined in (2.28) and (2.29) and according to 
(B.5), D^z; x, y) = D^z; x) - D^z; y), where 

D»(z;x) = -i9b [ d 2 i a e a ^d p arg(i-x) (4.21) 
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Inserting B^, and (4.15) in (4.20), we will get six terms: BPB, BPD and DPD, cor- 
responding to the [P^ + d^d"] part of (4.15) and BCB, BCD and DCD, corresponding 
to the e^ Qiy part of (4.15). One of these, namely DPD will be gauge dependent, as a 
consequence of the fact that the ^-dependent part of fig is not gauge invariant. As 
before, we will only consider the gauge independent part of this term. Let us evaluate 
now each one of the six terms above. 

The BPB term can be computed exactly as in (2.30), except for the fact that 
now the P°" A -proportional term between square brackets in (2.30) is replaced by 
^ r ~ 1 [ fc2 ^ fc 2 + g2) ] as one can infer from (4.15). We will arrive at the expression in the 
exponent of (2.33) but with (2.32) replaced by (4.16) in the square brackets in (2.33). 
Following the same steps that we took after (2.33) we arrive at 

r6 -0\x-y\ i 

\x-y\ 



BPB = \imnb 2 [ e l _ - - j-T ] + t / d 3 zd 3 z'A^(z; x, y)F(z - z')A^(z'; x, y) (4.22) 



where 

k 2 e -e\z~z'\ 

^-'>=^W 1 = - D sji=?i (4 - 23) 

Observe that the last term in (4.22) is surface dependent but no longer canceled by 
the renormalization counterterm S^LAJ. 

Let us obtain now the BCD term. This is given by( notice the factor two because 
there are actually two crossed terms) 



BCD = <\™/V p 6, fj (c) d 2 i„J T (o) dVa?W 9 (f- xjdfargtf- 

*F% (4-24) 
Using an identity similar to (2.31), namely, 

a e al3X e apx d p = e^e^ 1 {-U8 Xl + d x d~<) (4.25) 

we arrive at an expression identical to the exponent in (2.33) except for the prefac- 
tor and for the expression between square brackets which is replaced by jF~ 1 [ fc2( , fc j- +6) 2) ]- 
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Following the same steps we took after (2.33) we immediately obtain 

i _ p -0\x-y\ off 2 r 

BCD = 2n6 2 b 2 [ _ - - 1] + — J d 3 zd 3 z' A^z; x, y)G(z - z'jA^z'; x, y) 

(4.26) 

where G(z — z') = jF~ 1 [ fc2 ^_ g2 ] is given by the last equality in (4.16). Again, the last 
term in (4.26) is surface dependent. 

Let us now consider the DPD term. As we remarked before, it is gauge dependent. 
In order to get the gauge invariant part, let us take the <5 Mly -proportional part of (4.15). 
Then, using (4.21), we get an expression identical to (2.34), except for the prefactor 
and for the 5 3 (£ — rj) function which is replaced by ■7 r ~ 1 [ fc2 ^ 6)2 ] = G(z — z'), namely 

DPD = - j / d 3 zd 3 z'A^(z; x, y)G(z - z')A^(z'; x, y) (4.27) 

The gauge invariant part of DPD only contains the surface dependent term above. 
Let us evaluate now the crossed BPD term. This is given by 

2 

V 

p,8,m^0 ' J JT Xi {C) "" JT X .(C) 



BPD= lim -ieb 2 f2^j[ d 2 ^l d^ad^argd-x^argiff-Xj 



**S ^ J {^tl^^] - - ^1) <««) 

All the gauge dependent (derivative) terms vanish because F^ v a d G = 0. Let us 
observe now that 

F$ a e a/3x P$ = -a(d^e ua(S - d\^ aP ) (4.29) 

Inserting (4.29) in (4.28), we see that the second term vanishes because d 2 ^ / / d 2 rf . 
The first term gives 

BPD= lim - % Qb 2 Y.^A d 2 iJ d 2 r ]a dfarg{£-x l )dfarg{r)-x 3 ) 

P,S,m^0 =1 JT X .(C) JT Xj (C) 



where we used (4.16). 
Writing 



"^^[■j^it-T,) (4.30) 



Q2 

(4.31) 



k 2 + e 2 k 2 k 2 (k 2 + e 2 ) 
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and using the last equality in (2.32), we see that the first term of (4.31) when inserted 
in (4.30) leads to an integral identical to the one evaluated in Appendix D. Hence, 
using (D.2) and (4.31), we obtain 

2 

BPD = -inQb 2 €itj[arg(xi - £,-)] + K ( T ^ T y) (4.32a) 

or 

BPD = \imi2n6b 2 [arg{x - y) +arg{y-x) -2arg(e)} + K(T X , T y ) (4.326) 
where 

K(T x ,T y ) = lim iOW ]T ^ I d % I d 2 r) a d^arg(i- ^dfargtf- x 3 ) 

ij=1 JT Xi (C) JT Xj (C) 

xdfoF(t-T,) (4.33) 

where F(£ - rj) = F" 1 ^]^ }. 

Observe that as the previous terms, also BPD contains a surface dependent term, 
namely K(T x ,T y ), given by (4.33). 

Let us consider now the DCD term. Using (4.21) and the ^-dependent part of 
(4.15), we find immediately 

DCD = p lim^^ p e t e 3 jf {c) d% J t ^ d 2 Va d^arg(i - x^argtf- x s ) 

xe^ a e af3x e apx dpF(£-r)) (4.34) 
where F(£ - rj) = J r ~ 1 [ fc2(fc l +g3) ]. Observing that 

e ^a e afl\ e ap\ d (t) = e <*P»^ _ e ^^ } (4 35) 

we see that the second term gives a vanishing contribution because d 2 £ M / /d 2 £ a . The 
first term is identical to (4.33) and we get 

DCD = ~K(T X ,T V ) (4.36) 
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Let us obtain, finally, the BCB term. This is given by 

i6b 2 2 r f -> 

BCB = lim^ — J2 e i e i J T {c) d % J T {c) d 2 Vadi°arg{£- x^df arg{rf - Xj) 

xFfcrFtf^djPFit-rj) (4.37) 

where again F(£ - rj) = ^[ p^y ]- 
We now have the identity 

-e^d^dj^ (4.38) 

Again, the contribution of the 2 nd term vanishes because d 2 ^ / /d 2 r] a . In Appendix 
E we show that the contribution of the 3 rd and 4 th terms also vanish. We also show 
in Appendix E that the 1 st term gives 

BCB = ~\ BPD (4.39a) 

or 

BCB = \im-iix6b 2 [arg(x - y) + arg{y- x) - 2arg(e)] - ^K{T x ,T y ) (4.396) 

Collecting the six terms given by (4.22), (4.26), (4.27), (4.32), (4.36) and (4.39), we 
see that the surface dependent part of the last three terms, namely K(T x ,T y ) exactly 
cancel! On the the other hand, we can see, by using the fact that 

^WT^ + ^FT# ] = = 63{z ~ z>) (440) 

that the sum of the surface dependent pieces of the first three terms precisely add up 
to Sfll/LJ and therefore are exactly canceled by the renormalization counterterm in 
(4.20)! Just the surface independent part of the above six terms contributes to the 
correlation function and we get: 

e -6\x-y\ i _ e -0\x-y\ 

< fie, R (x)fil R (y) >= exp{vro 2 [^^] + 2n6 2 b 2 [ _ ] 

\x y\ \x y\ 
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+i7c9b 2 [arg(x — y) + arg(y — x)]} (4-41) 
where we defined the renormalized fig as 

nb 2 

He,R = lim fig exp[— + nb 2 9 2 - in9b 2 arg(e)} (4.42) 

Observe that for 9 — > (4.41) reduces to (2.37). Eq. (4.41) is multivalued up to a 
factor exp[i27r(7r#6 2 )], indicating the fig has statistics = -ir9b 2 . This is in agreement 
with the previous result (4.9). 

We see that fig is indeed a local field. Notice that the correlation function of fi 
would be given by exp[BPB + BCB], indicating therefore that fi would be a nonlocal 
(surface dependent) field in the MCS theory. The same would be true for £' (even in 
Maxwell theory) whose correlation function would be given by exp[DPD + DCD]. We 
see that the operators studied in [24] are nonlocal in the MCS theory and therefore 
cannot be used as bona fide dual operators. 

The above study reveals how delicate and stringent are the conditions necessary 
for the obtainment of local dual operators. 

Observe that 

< fJ>8R{x)v*8 iR {y) > |:E constant ^ 

indicating that fig does not create genuine vortex excitations in the MCS theory. In 
section 5 we will see how this result is changed when we add a Higgs potential to the 
MCS theory. 

4.5) Mixed and Composite Anyon Correlation Functions 

Let us study here the mixed a — fi</> correlation function, from which we will obtain 
the composite anyon tp correlation function. 
Combining (4.13) with (4.19), we can write . 

< a(x 1 )fig(x 2 )a*(y 1 )fi* e (y 2 ) >= lim Z' 1 f DW„exp{- / d 3 z[^-W fl [P^+9C^+G^}W l/ 
+Cp(z;x 1 ,y 1 )W + W,[P^ + 9C^]A l/ (z;x 2 ,y 2 ) + ^(z 2 ,y 2 )]} (4.43) 
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Integrating over W^, we get 

< (T(x 1 )ne{x 2 )a*{y 1 )if e (y 2 ) > - lim exp{- f d 3 zd 3 z'[C^z;x 1 ,y 1 ) + B^z; x 2 , y 2 ) 

+D^(z; x 2 , y 2 )] [C u (z; x 1 , yi) + B u (z'; x 2 , y 2 ) + D u (z'; x 2 , y 2 )\ 

x[P^ + BC" + G» v Y\z - z') - S R [A,(x 2 ,y 2 )]} (AAA) 

where was defined in (2.39-40), in (2.28-29) and in (4.21). Inserting expres- 
sion (4.15) for \P^ V -\-0C' IU + G fiu ]~ 1 , we will get, in addition to the six terms computed 
in (4.4) six more terms involving C M , namely CPC, CPB, and CPD, corresponding 
to the [P^ + d»d u ] part of (4.15) and CCC, CCB and CCD, corresponding to the 
e mv part of (4.15). As before, we only consider the gauge invariant part of eventually 
gauge dependent terms. 

The CPC and CCC terms were computed in section 4.3. Their sum is given by 
the exponent in (4.17b). The CPB term was computed in section 2.5. It is given by 
(2.51). The gauge invariant part of the CPD term is given by 

CPD = lim ab9 X; Utj I d\ f d 2 Va [e ik d®arg(£ - n) + % In \£ - f t \] 

Me-0 JRr x (C) JT Sj (C) 

. e -e[\t-r,\ 2 +\A 2 ] 1/2 

where we used (4.16), and introduced the ultraviolet regulator e. ( We use the same 
convention for r, and Sj as in section 2.5). Using now (2.43a) and the results of 
Appendix C, we find 

CPD = ling -f ± e M - eM ~ elir '-f + lCm ] (4.46a) 

or 

CPD = -nab[e- elxi - y21 + e ~ elx2 - yi1 - e^ 1 ^ 21 - e^ 1 ^ 21 ] (4.466) 
Let us consider now the CCD term. This is given by 

CCD = lim iab6 2 £ [ d 2 i [ d 2 Va [e ik d^arg(£- f z ) + dL In \£- f,|] 

P^O . j=l jRr^C) Jt Sj (c) 
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xdfarg(ff- Sj)e a ^e ipx d^ F(f - rj) (4.47) 

where F(£ — rj) — jF~ 1 [ A;2( - fc2 1 +6)2 ^ ]. In Appendix F, we show that CCD = 0. 
Let us finally evaluate the CCB term, which is given by 

CCB = lim -aW ]T J r ^ cP^d^arg^- n) + 3fo In |f- f,|] 

x^argtf- sj)F^ x j*d®F(Z - V ) (4.48) 
where F(£ — 77) is the same as in (4.47). Observe now that 

Iff) xe ipX df = e*%% - ^ Q <)% ( 449 ) 

As is shown in Appendix C, the only possibility for the ^-integral in (4.48) to be 
nonzero is the expression between brackets (in (4.48)) being contracted with a deriva- 
tive d l . We could try to force the appearance of such a derivative by using an identity 
such as (2.43). Such identity, however, would also be contracted with (4.49) and we 
immediately see that the derivative part of the identity would vanish, because of the 
contraction with the e's and d p . We conclude, therefore, that CCB = 0. 

Collecting all terms contributing to the mixed correlation function, we get 

< vMtoAxtWyJtoAv*) >= exp{^V e ^-^ - 1]+ 

p -8\x2-y2\ 1 _ p ~9\x2-y2\ 

nb 2 [j r ] +2tt0 2 & 2 F— 1 -Trftre-^i-wl + e -% a -wl _ e -»l*i-*2| _ e -%i-wli 

F2-2/2I \x 2 -y 2 \ 

+inab[arg(x 1 - x 2 ) + arg(y x - y 2 ) - arg{x x - y 2 ) - arg(yi - x 2 )\ 

+m9b 2 [arg(x 2 - #2) + arg(y 2 - x 2 )]} (4.50) 

Observe that this expression reduces to (2.52) when 9 — > 0. The correlation 
function for the composite anyon field (2.53) can now be easily obtained by taking 
the limit x\ — ■> rc 2 and yi — > y 2 in (4.50). Introducing the same renormalization factor 
as in (2.53), we get 

< <p(x)<p*(y) >= exp{[— - 27ra6][e- 9 l a; -^ - 1] + tt6 2 [^ -] 
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+27rfl 2 6 2 [ 1 6 _ X V ]-i7r{ab-6b 2 )[arg(x-y) +arg{y-x)]} (4.51) 
\x y\ 

Notice that the multivaluedness of < if if* > corresponds to a phase exp[i2n(Trab — 
ir9b 2 )]. This is in agreement with (4.11) and indicates that if has statistics S v = 
nb\a — 9b\. Observe that for 9 — > 0, (4.51) reduces to (2.54), which we obtained in 
Maxwell theory. 

5) The Maxwell-Chern-Simons-Higgs Theory 

5.1) The Unbroken Phase 

Let us study now the Maxwell-Chern-Simons-Higgs (MCSH) theory which is es- 
sentially the Abelian Higgs Model plus a Chern-Simons term: 

£ = ~\Kv + l^WpdaWp + \D,<j>\ 2 - g 2 <t>*<t> - ^(0» 2 (5.1) 

As in the AHM, we have two phases, according to whether: i) g 2 > 0, < >= or 
ii) g 2 < 0, < >= k ^ 0. Shifting the Higgs field around n (for g 2 < 0), as before, 
will generate a mass term (M = en) for W^. 

Let us introduce here the operators a and fig considered before. The equal time 
commutation rules of the MCSH theory are the same as in (4.5) and therefore all the 
commutators evaluated in section 4.2 remain valid here. 

The same remarks concerning the surface independence of \x in the AHM may be 
applied to fig here and therefore we conclude that (2.7) is the appropriate surface 
renormalization counterterm for fig in the MCSH theory too. 

Let us consider firstly the symmetric phase (g 2 > 0, < >= 0). The /i^-correlation 
function is given by 

< fie(x)fi* e (y) > lim Z' 1 f DW A1 exp{- f d 3 z[ V M [P^ + BC^ + G^]W V + |£>^| 2 + 

V{4>) + WyP^ + QC^\A v {z- x, y) + ±A%,(z; , x, y)]} (5.2) 

As in the AHM, < figfi* e > will be given by (3.3). Again, only two-leg graphs will 
contribute to the large distance limit in (3.3) [17]. Also here, we are going to make a 
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loop expansion for the computation of A(x, y) in (3.3). As in the AHM, the leading 
contribution will be given by the graph of Fig. 4b. The only difference now is that 
the W^-propagator is given by (4.15) and the external field is B^(z; x, y) + D^z; x, y), 
where is given by (2.28) and by (4.21). We immediately see that the graph of 
Fig. 4b, for these external fields and propagator is identical to the first term in the 
exponent in (4.20). It follows that 

lim < HeA^^eAv) >mcsh=< toA x )P*6,nhi) >mcJ X constant ^ 

\x-y\— >oo 

(5.3) 

where < ^e,R^e,R > mcs is given by (4.41) and we renormalized /ig as in the MCS 
theory. 

We conclude that also in the ordered phase of the MCSH theory the fig operator 
does not create genuine vortex excitations. We will see that this is no longer true in 
the broken phase of the theory. 

Let us evaluate now the order correlation function < a a* >. This is given now by 

< a(x)a*(y) > lim Z' 1 f DW„exp{- / d 3 z[-W^[P^ + 6C^ + G> XV \W V + \D^\ 2 

+V(</>) + C lt (z;x,y)W' i ]} (5.4) 

where Cy^zjx, y) was defined in eqs. (2.39) and (2.40). As in the AHM, < aa* > is 
given by (3.6). Making exactly the same approximation we did in section 3.2 we see 
that r(x, y) will be given by the graph in Fig. 4b but with C M as the external field 
and (4.15) as the H^-propagator. It is easy to see that this is given by the exponent 
in (4.14). We readily conclude that in lowest order, the long distance behavior of the 
gauge invariant part of < aa* > in the MCSH theory is given by (4.17), namely, 

lim < a R (x)a* R (y) > MCSH =< a R (x)a* R (y) > M cs X 1 (5.5) 

\x-y\-*oo 

where < a R a* R >mcs is given by (4.17c) and we renormalized a as in the MCS theory. 

We now see that the behavior of < cr R a* R > indicates charge screening in spite of 
the fact that we are in the symmetric phase. This naturally happens, because as in 
the MCS theory, the W^- field acquired a mass through the Chern-Simons mechanism. 
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As before, the long distance behavior of the composite anyon field could be ob- 
tained by just exchanging B^z; x, y)+D^(z; x, y) by B^z; x, y)+C fl (z; x, y)+D fl (z; x, y) 
in (5.2). We would obtain (4.51) for the long distance behavior of < ipip* >. 

5.2) The Broken Phase 

Let us investigate now the broken phase, where < >^ 0. Shifting the Higgs field 
around the vacuum value, as before, we can see that a mass term will be generated for 
Wp. The quadractic lagrangian, however, contains a Chern-Simons term in addition 
to (3.8). The euclidean W^- propagator will be now given by 

D (k) P " U{1 + ^ ~ 6e ^ ka I W (56) 
^ ] (k 2 + M 2 ) 2 + k 2 6 2 k 2 {£k 2 + M 2 ) y ' 

Observe that this reduces to (3.9a) for 9 — > and to (4.15a) for M — > 0. 

Let us evaluate the order correlation function < aa* >. This is given by (5.4) 
and (3.6). Making the same approximation as before, the leading contribution to the 
long distance behavior of < aa* > will be given by the graph of Fig. 4b with C M 
as the external field which, as we pointed out before, leads to an expression similar 
to (4.14). The only difference is that now we must use the M ^ and 9 ^ 
propagator (5.6) instead of (4.15). As in the previous cases we are going to compute 
the gauge independent part of < aa* >, namely the 5^ and & av proportional terms. 
As in the MCS theory, the e^^-proportional term vanishes (see Appendix C). The 
<5 MI/ -proportional contribution, namely 

5 ^-ir k2 + M " i = r ikk ^M r k2 + M2 1 (5 7) 
[ {k 2 + M 2 ) 2 + k 2 9 2i 2n 2 Jo \x\ l {k 2 + M 2 ) 2 + k 2 9 21 y ' 

can be calculated with the help of the identity 

sin A; la; I „ ,, r cosfclxl, x^x v r sinA;|x| , ,. M nN 

$ — n — = dd [ r^- 1 +-r-nr — —^-kcosklxl] (5.8) 

\x\ k \x\ 2 \x\ 

which generalizes (2.43) when inserted in (5.7). Proceeding exactly as we did to go 
from (2.42) to (2.44), we see that the contribution from the last two terms in (5.8) 
vanishes. The contribution from the first term gives 

< a(x)a*(y) ^J^ 00 e xp{47rV[/(:r - y) - 1(0)}} (5.9a) 
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or 



< a R (x)a R (y) > l ^°° exp{An 2 a 2 I(x - y)} (5.96) 



where 



' f°° k + M 

I{x - y) = 2^Jo dkcosklx - yl[ WTM^Tk^ ] (5 - 10) 

and a R = (jexp[27r 2 a 2 /(0)]. Observe that 

-M\x-y\ 



liml(x-y) = 6 - — (5.11) 

and (5.9) reduces to (3.11) in the AHM. 

It follows from the Riemann-Lebesgue lemma that lim.| :r _ y |_ i . 00 I(x — y) = 0. Hence, 
< a R (x)a* R (y) > — > 1 . This result indicates that also in the broken phase of the 
MCSH theory the charge is screened. Actually it is doubly screened, both by M and 
by 9, namely by the Higgs mechanism and by the Chern-Simons term. Notice that 
for M — > and 9 — > 0, we have 

lim I(x-y) = — [— - \x-y\] (5.12) 

and (5.9) reduces to the expression (3.7) found in the symmetric phase of the AHM, 
where charge screening no longer occurs. 

Let us consider now the disorder correlation function < Ho^o >■ According to 
(5.2), we can write < > exactly as in (3.12), with the only difference that now 
SlW^, 0, D^4>] is the action corresponding to (5.1). Making, as before, the change 
of variable — > + A^(z; x;y), we get an expression identical to (3.13), where 
again, S is associated with (5.1). The correlation function < fie^e > therefore, will 
still be given by (3.14), the only difference being that the W^- propagator used in the 
computation of Feynman graphs must be (5.6) instead of (3.9a). As before, making 
an expansion in loops and in powers of M 2 , it follows that the leading contribution 
to A(x,y) is given by the graphs of Fig. 6. Using (5.6), it is easy to see that the 
sum of the gauge dependent graphs of Fig. 6a identically vanishes for all values of 
the gauge parameter £. The first contribution (of order 0(^)) is given, as in the 
AHM, by the graphs of Fig. 6b. These are independent of 9 and we see that in lowest 
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order the /^^-correlation function behaves asymptotically at large distances exactly 
as (3.18). The ^-dependence will only be introduced by 1-loop corrections, through 
the W^-propagator. 

The result < fj,e jR (x)(j,0 R (y) > — > in the broken phase of the MCSH theory 
shows that ho,r creates true topological charge bearing excitations in this case. The 
mass of these quantum solitons is M v = nM 2 b 2 , within our approximation. 

Let us remark finally that the (charge and magnetic flux bearing) composite anyon 
operator (p would be the creation operator of the quantum excitations corresponding 
to the classical electrically charged vortices of ref. [25]. 

6) Conclusions and Remarks 

The order-disorder duality structure was exploited in order to introduce dual op- 
erators corresponding to a local U(l) symmetry and carrying respectively, charge and 
magnetic flux, in some 2 + 1 dimensional field theories involving a vector field. The 
conditions for the obtainment of local order (a), disorder (vortex) (/i) or anyon (ip) 
operators must be carefully analyzed in each case. Local, surface independent cor- 
relation functions of these operators can be obtained only after some very stringent 
requirements are met. 

In the absence of a Chern-Simons term, namely, in the Maxwell and in the Abelian- 
Higgs theories, the behavior of < a > and < fi > which can be infered from the long 
distance behavior of the a and fi correlation functions is analogous to the behavior of 
the corresponding operators in the Ising model. The Maxwell and symmetric Higgs 
phases correspond to the disordered phase of the latter and the broken Higgs phase 
corresponds to the ordered phase of it. The ordered phase contains vortex excitations 
whose mass is explicitly evaluated. The disordered phase may be viewed a vortex 
condensate. Massive anyon excitations occur in all cases. The screening of charge in 
the broken phase is clearly exposed by the behavior of the a-correlation function. 

In the presence of a Chern-Simons term, namely, in the Maxwell- Chern-Simons 
and Maxwell-Chern-Simons-Higgs theories, we have two concurrent mechanisms of 
charge screening (mass generation for the vector field). The behavior of a-correlation 
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functions expresses the Chern-Simons induced charge screening in the MCS and sym- 
metric MCSH phases. In the broken MCSH phase we can infer from the long distance 
behavior of < aa* > the action of both mechanisms of charge screening, namely, the 
Chern-Simons and Higgs mechanisms 

Massive vortex excitations are still only seen to occur in the broken phase of the 
MCSH theory. We therefore can clearly distinguish the process which generates a 
mass for the vector field from that which generates a mass for the vortices. The first 
one can occur either via the Chern-Simons or Higgs mechanisms, while the second 
one is only induced by the latter. 

In the presence of a Chern-Simons term we show that the statistics is no longer 
determined by the product (charge x flux) but rather, by where Q = Q — 6&. 
As a consequence, the exclusively flux bearing disorder (vortex) operator \i$ is seen 
to be anyonic. 

The charge and magnetic flux carrying quantum vortex excitations corresponding 
to the classic solutions of the MCSH theory found in [25] will be described by the 
composite anyon operator if = cr/^. 

Let us remark that even in the presence of a Chern-Simons term, an Ising-like 
behavior for the order and disorder variables could be obtained, by introducing a new 
order operator a' = <7<f), where <fi is the Higgs field. This could be done naturally in the 
Abelian Higgs and MCSH theories but could also be easily achieved in the Maxwell 
and MCS theories by considering the weak coupling limit of them with a complex 
scalar field (without self-interaction, for instance). The behavior of < a' > in the 
MCS and MCSH symmetric and broken phases would be the same as that of < a > 
in the Maxwell and Abelian Higgs (symmetric and broken) phases. Furthermore, the 
short distance behavior of < a' a'* > in the Maxwell and MCS theories would have 
the singularity one would expect usually. The main advantage of the operator a is 
that it provides a reliable order parameter which works well whether or not a Higgs 
field is present in the theory. 
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Appendix A 

Let us show here the Cauchy-Riemann equations involving the real and imaginary 
parts of the analytic function ln(z — x) that we use in this work. 

The first one is (we will be in euclidean space throughout this Appendix) 

POO 

d ia rg{z- x) = -e ij dj In \z - x\ + 2n / d^e ij 5 2 (z - f ) , i = 1, 2 (Al) 

The line integral corresponds to the singularity of the derivative of arg(z — x) 
along its cut L, chosen to be the straight line going from x to oo along the z 1 axis. 
From (A.l) we get a second equation by contracting with e u : 

roc 

e ij diarg{z- x) = -d j In \z - x \ + 2vr / d£ ] 5 2 {z - f) , i = 1, 2 (A2) 

Jx,L 

From (A. 2) we obtain 

e ij didjarg{z- x) = 2tt5 2 (z -x) + 2n f°° dCdf ] 5 2 {z - f ) (A3) 

Jx,L 

where we used the fact that did 1 In \ z — x\ = 2n5 2 {z — x). 

Appendix B 

Let us derive here eqs. (2.13) and (4.6). We start with (2.6) and observe that 
= d^A u - dyAp = / d 2 i u arg(£- x)d^{z - £) - (// <-> v) 

+ / d^e^argii- x)5\z - (B.l) 
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In this equation, the second term comes from the discontinuity of at C(x), the 
boundary of T X (C) (see Fig. 1). Inserting (B.l) in (2.6), we get 

fi(x; C) = exp{tb f d^d^arg^-x) -% I d^e^W^ 'arg(£-x)} (B.2) 

JT X (C) I JC(x) 

Observing that d fl W tlu = e fiap e uXp d a d\W fl and integrating by parts the first term in 
(B.2), with the help of Stokes theorem we see that the boundary term exactly cancels 
the last term in (B.2). The remaining term gives 

fi(x; C) = exp{-ib [ d 2 £ u W u d,arg(£- x)} (5.3) 
Jt x (c) 

Noting that d 2 C, u only has the zeroth component and remembering that E l = W l °, eq. 
(2.13) immediately follows. 

In order to obtain eq. (4.6) for fig, let us consider £', eq. (2.10b). Inserting (B.l) 
in (2.10b) (we can write this equation in terms of A pu divided by two) we get 

E'(x; C) = exp{ia f d 2 ^ a/3 d a Wparg(£ - x)+iai d^W"arg(£- x)} (BA) 

Integrating by parts the first term with the help of Stokes theorem we again find 
that the boundary term exactly cancels the last term. The remaining term yields the 
following expression: 

T,'(x; C) = exp{-ia f d 2 ^ aP W p d a arg(£ - x)} (B.5) 

This equation, for a = 9b together with (B.3) leads to expression (4.6) for fi e . 

Appendix C 

Let us demonstrate here the following two useful results: 

i?! = lim / d 2 i[e ij diarg(i- x) + & In \£ - x\]djA(x°, |) = 2tvA(x°, x) (C.l) 
and 

# 2 = lim / d 2 £e ij d i arg{£-x)d j K{x°,x) = 2itA(x°,x) (C.2) 

p,S-*0j Tx (C) 

for an arbitrary function A(x). 
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Let us start with (C.l). We may write Ri as (we ommited the argument x°) 
R-y = lim / d^d^argit - £)A(f)] - e^d^arg^- £)A(f) 

+di[di In |f - £ |A(f)] - 9,9, In |f - x |A(|) ] (C.3) 
Using (A. 2) or (A. 3), we see that the sum of the 2 nd and 4* h terms is: 

2 nd + 4 th = lim 2tt / d 2 £ H di 1 j df ) 5 2 {i-r 1 )A{£) (CA) 

P,<5^0 JR X (C) Jx,L 

= lim 2tt / dmd\ v) A(ff) = lim 2tt[A(£ + pf ) - A(f)] = (C.5) 

p,5^0 Jx,Lr\R x (C) P^O 

where the last integral is performed over the section of L contained in R X (C). Observe 
that (A. 3) is not zero in a cut plane like R X (C). 

It would also be different from zero in a punctured plane, for instance (this is the 
case considered in [11]). 

The I s * and 3 rd terms in (C.3) may be evaluated with the help of the Stokes and 
(two dimensional) Gauss theorems, respectively. The result is 

R x = lim {- / dC&argd- x) + e^d 3 In |f - x |]A(f)} (C.6) 

In this expression, C p is the arc of circumference of radius p and Cs, the two straight 
lines, in Fig. 3. is the arc of circumference with infinite radius closing C at infinity. 
We immediately see that the contribution from Cs is zero because the integrand is 
orthogonal to the integration element along Cs- The contribution of C^ also vanishes 
because of (A.l). The only contribution to (C.6) comes from C p . Inserting (A.l) in 
(C.6) we get , 

i?! = lim 2vrA(x°, x + pf ) = 27rA(:r , x) (C.l) 

thus establishing (C.l). 

Let us observe that replacing djA by a vector Aj (not a derivative!) in (C.l) 
would yield a zero result. Indeed, inserting (A. 2) in (C.l) we would obtain an integral 
identical to (C.5), with djA replaced by Aj. This would vanish at p — > for any 
regular Aj. The boundary terms, eq. (C.6) which are only present when Aj = <9jA 
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are responsible for making R\ ^ 0. This is the reason why expressions of the type we 
found in the last terms of (2.43) and (5.8) do not contribute to correlation functions. 

. . — * 

Let us turn now to (C.2). Observe that e^djdiarg^ — x) = over the surface 
T X (C)), Fig. 1. Using this fact, we can make the dj derivative total in (C.2) and then 
use Stokes theorem to get 



R 2 = lim & dCdiarg(i-x) (C.8) 

p,5->0 JCoo+C s +C p 

Here C p is the arc of circumference of radius p and C$, the straight lines in Fig. 1. C^ 
is the infinite radius arc of circumference closing C at infinity. The contribution from 
Cs again vanishes because the integrand is orthogonal to d^ 1 along it. The contribution 

— * 

from Coo vanishes for A(£)-functions going to zero at infinity. This is true for all cases 
considered in this work. Even for the function between square brackets in (2.43b) 
this is seen to be valid by adopting the following procedure before taking the limit 
x — > oo: we put the system inside a circular box of radius R = ^. We immediately 
see that the above mentioned function will vanish on the boundary of this box. Then 
we make R — > oo as x — > oo. 

The only nonvanishing contribution to (C.8) therefore comes from C p . Making a 
Taylor expansion in A(£) around £ = x, we see that in the limit p — * only the zeroth 
order term is nonvanishing and we obtain: 

R 2 = lim A(x°,x) [ d?diarg(£ - x) = 2nA(x°,x) (C.9) 

P,<5^0 Jc p 

This establishes (C.2). 

In eq. (4.7b), we have an expression similar to (C.2), the difference being that 
the boundary of the integration region now contains two curves: C x and C y . Using 
the Stokes theorem with dj as the total derivative in C x and with in C y , eq. (4.8) 
immediately follows from (C.2). 

Appendix D 

Let us demonstrate here the following two useful results: 

R 3 = lim / d^le^d^argiff-^ + d^hilff-yl] 

px,Sx,Py,Sy^O JRy(C) 
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x J Tx(c) d%dParg(i- ^[-^—fl = narg(y- x) (DA) 



and 



= 7r[arg(|/ — rr) + arg(x — y)] (-D-2) 
Let us first demonstrate the following result that we will need later 

It may be written as 

h = lim / dpp / rfrr J ff[L 2 + r 2 -2Drcos(v?-v?o)]~ 3/2 (DA) 

Px,S x ^0Js x Jp x 

where L = \rj — x\, ipo = argifj — x), H = L cos# and D = L sin6>, where 9 is the 
angle the 3-vector r\ — x makes with the £ 3 -axis in 3D-euclidean space. 
Evaluating the r-integral [26], we obtain 

r2n-S x H 

L = lim / dip ip- — ; r (D.5) 

Js x L - D C0S(V9 - (fo) 

Performing now the (^-integral [27], we get the last equality in (D.3). 

Let us consider now the last integral in (D.l) and (D.2) (observe that v = % — 1,2 
in (D.2), because d 2 r] a has only the 3-component different from zero. Also fi = 3 in 
(D.l) and (D.2) for similar reason): 

Integrating by parts and using the 2_D-Gauss theorem, we may write (D.6) as 

Using (D.3) in the first term of (D.7), inserting the result in (D.l) and (D.2) and 
using (C.l) and (C.2), respectively, we see that the first term of (D.7) contributes a 
term narg(y — x) to (D.l)and (D.2). 
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We will now show that the second term (boundary term) in (D.7) only contributes 
to (D.2). Indeed, using (A. 2) we see that in the limit 5 X — > 0, the first expression 
between brackets in (D.l) is orthogonal to e jl d£ j along the two straight lines in C(x). 
The contribution of the arc of circumference with radius p vanishes for p — > 0. Hence, 
only the first term of (D.7) contributes to (D.l) and the result of the previous para- 
graph establishes (D.l). 

Let us compute now the contribution from the second term in (D.7) to (D.2). 
Inserting this term in (D.2), we obtain the boundary contribution to R^: 

Observe that the second integral again has an expression like (D.6). Inserting (D.7) 
in (D.8) we see that the second term vanishes because d£ j is orthogonal to e kj drj k 
along the straight lines of C x and C y . As before, the contribution from the arcs of 
circumference vanish for p x ,p y — > 0. Using (D.3) we see that the first term of (D.7) 
gives the following contribution to (D.8): 

R? =-7; Km „ / dt 3 arg(£- x)df ] arg(£- y) 

= - Inn A d^[df ] arg(i- x)]arg(£- y) = irarg(x - y) (D.9) 

where we used (C.9). Combining this term with the contribution from the first term 
of (D.7) to i? 4 , we immediately establish the last equality in (D.2). 

Let us finally mention that actually a more general result could be established for 
(D.l) and (D.2). Before taking the limit p x ,5 x — > in (D.l) one would obtain R 3 = 
\arg{y — x)Q(y,T x (C)) , where Q(y,T x Q) is the solid angle determined by the point y 
and the surface T X (C). For p x , 5 X — > , fl(y, T X {C)) — > 2n and we recover (D.l). Also, 
without taking the limits in (D.2) one would obtain i? 4 = ^[arg(y — x)Q(y,T x (C)) + 
arg(x — y)Q(x,T y (C))], which reduces to (D.2) for p x ,5 x — > and p y ,S y — > 0. These 
results would be relevant for the computation of correlation functions of the loop 
dependent operators p(x;C) (finite p,S). 
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Appendix E 

Let us demonstrate here the result (4.39) for the BCB term which is given by eq. 
(4.37). 

We start showing that the contribution from the last two terms in (4.38) vanishes. 
Insertion of the 3 rd term of (4.38) in (4.37) yields an integral 

H(xi,Xj) = lim / d 2 ^ d 2 rj a d^arg(£ - x^d^ arg(ff - Xj) 

p,<5->0 JT Xi (C) JT Xj (C) 

xe^d^d^-rj) (E.l) 

Making the change of variable £ <-> rj and exchanging a <->//, f3 <-> v, and using the 
fact that dr£\F(£ — rj) — —d/^F^ — rj), we immediately conclude that the 4 th term 
in (4.38) would yield an integral H(xj,Xi). Hence, summing over % and j in (4.37) we 
see that the contributions from the last two terms in (4.38) are identical. Let us now 
show that these contributions actually vanish. Using the fact that d 2 ^ and d 2 r] a only 
have the 3-component different from zero and noting that arg(ff— Xj) — along the 
surface T Xj (C), we may write 

H( Xi , x 3 ) = lim / d 2 Ze kl d®arg(£- x,) 

p,S^0 JT Xi (C) 

xd (O d hf , fr&Wargtf-fyFit-Tj)] (E.2) 

JTxj (C) 

Using now the Stokes and 2D-Gauss theorems, we get 

H(x i ,x j )= lim / d^d^argg-^l drfe ji d^arg(ff - Xj^Ffc - 77)] 

(E.3a) 

= lim o 2vr£ ( ^drfd^ In \ff - Xj\[df Xi) F(xi - rj)] (EM) 

where we used (C.9) and (A. 2). 

The contribution from the arc of circumference of radius p (Fig. 1) to (E.3b) 
vanishes because drf is orthogonal to the integrand along this curve. The contributions 
from each one of the straight lines in C(x 3 ) (Fig.l) cancel each other in the limit 
S Xj — > because the integrand is regular along the cut. Hence we conclude that 
H(xi,Xj) = 0. 
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The only nonvanishing contribution to (4.37) comes from the first term in (4.38). 
Insertion of this term in (4.37) yields the following integral 

J(xi,Xj)= lim / d 2 L / d 2 r} a d^arg(£- Xi)d^arg(ff- xj) 

p,5->0 jT Xi (C) JT Xj {C) 

where F(f - rj) = ^[ ^J^ ]- 

Before evaluating (E.4), let us demonstrate the following results that we will need 
later: 

/ d ^d,arg(i-x)^p^f^-y\) = (E.5a) 
j t x {c) |4 - y\ 

and 

/ d 2 i d^argd- £)^^/(|£ - y\) = (E.Sb) 
Jt x {c) |4 - y\ 

where / is an arbitrary function. 
We may write (E.5a) as 

J P J s r r [ [D 2 + r 2 -2Drcos(^-^ )] 1/2 

(£.6) 

where L, D and </? were defined in Appendix D. The angular integration in (E.6) is 
easily seen to vanish, because the contribution from the first two quadrants exactly 
cancels the one coming from the last two. This establishes (E.5a). Eq. (E.5b) 
immediately follows if we note that it will be given by an expression like (E.6) but 
with | replaced by — 

Let us now return to (E.4). Observe that d?s only has the 3-component different 
from zero (because it is contracted with d 2 r] a ). Since arg(ff— Xj) does not depend on 
if, we may make a total derivative and use Gauss theorem on the closed surface S, 
consisting of T x . (C) UT^C) UT Xj (C). Here T^C) is a copy of T x . (C) but at r] 3 = oo 
and T Xj (C) is the surface paralell to the r/ 3 -axis, connecting T Xj (C) and T^C), which 
is obtained by translating the curve C(xj) from rf = x 3 to rf = oo along the ?7 3 -axis. 
For p, 5 — > of course, T Xj {C) collapses to zero. After using the Gauss theorem we 
may write (E.4) as, 

J(xi, Xj )= lim / d%[ d'rje^d^argd-x^dfargifi-Xj) 

P,6^0JT Xi (C) JV(S) 
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xd^i-D.F^-rj)} (E.7) 

where V(S) is the volume bounded by S (note that d 2 rj a points inwards V(S)). To 
arrive at (E.7), we used (E.5b) and the fact that Oarg(rf — Xj) = on the volume 
V(S). Observe that the contribution from T 0O (C) vanishes because of the boundary 
conditions at infinity we impose on the function F (see Appendix C). One may also 
show that the contribution from T Xj (C) vanishes for p, 5 — - > as should be expected. 

Observe now that dj^ also only has the 3-component different from zero in (E.7). 
This is so because it is contracted with e vlf5 and v,f3 = 1,2 since these indexes are 
contracted with derivatives of arg functions. We therefore can pull <9 7 as a total 
derivative in (E.7) (after using the fact that dj^ = ~dj^) in the same way we did 
with &fo in (E.4) and use Gauss theorem backwards to obtain 

J(xi,Xj) = lim / d 2 ^ d 2 r} a d^arg(£- x^dW arg(rj - Xj) 

x<r#dfa[-nF(t- v )] (E.8) 

Again we used the fact that the contribution from (C) is zero for the same reason 
above. The contribution from T x . (C) also vanishes because the integrand in (E.8) is 
orthogonal to d 2 r] a along T Xj (C). 

The integral Jfaxj) in (E.8) is identical to the one in (4.30) (-OF = J 7-1 [(A; 2 + 
d 2 )^ 1 ]), therefore, insertion of (E.8) in (4.37) immediately demonstrates (4.39). 

Appendix F 

Let us demonstrate here that (4.47) is equal to zero. Using the identity e al3X e tpx = 
5 m 5 f3p — 8 ap 8^ 1 , we see that the first term vanishes because d 2 rf is orthogonal to the 
expression between square brackets in (4.47). The second term, when inserted in 
(4.47) yields the integral 

L( Xi , Xj) = lim (-) / d 2 i [ d 2 Va [e ik d^arg(i- f.) + Sf Q In |f - n\] 

p,6^0 JR ri (C) JT Bj (C) y ' 

xd^argifj-s^Fit-r)) (F.l) 



53 



Observing that &fo = —d^ only has the 3-component different from zero and using 
the Gauss theorem exactly in the same way we did in Appendix E, we may write (F.l) 

as 

L(x u Xj ) = lim / d 2 i f d 3 V [e tk d ( farg(£- n) + flfo In ||*- n\] 

P,o^0JR ri (C) JV(S) 

xd^arg(ff- ^[-D,^ - rj)} (F.2) 

where V(S) is the volume introduced in Appendix E. 

As we have shown in Appendix C, the only possibility for the ^-integral to give a 
nonzero result is the ^-integral in (F.2) producing a total derivative. Using the identity 
given by (2.43a), for [— OF] in (F.2), we see that the ^-integral of the derivative term 
vanishes according to (E.5a). Hence we conclude that lim Pr . )( 5 r .^ L( Xi , Xj) = and 
therefore CCD = 0, in eq. (A.47). 

It is worth mentioning that for the special case F = J r-1 [p], —OF = 5 3 (£ — rf) in 
(F.2) and then, we would obtain L(xi, Xj) = 2narg(fi — Sj). This, however is not the 
case of the theories studied here. 
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Figure Captions 



Fig. 1: Surface used in the definition of fi(x; C). 
Fig. 2: Surface used in the definition of a(x). 

Fig. 3: Surfaces used in the definition of the external field A^x; C). 
Fig. 4a: Vertex involving the external fields _B M , C M or D^. 

Fig. 4b: Leading graph contributing to the long distance behavior of correlation 
functions in the AHM and MCSH theory. 

Fig. 5: Vertices relevant for the evaluation of < fifx* > in the broken phase of the 

AHM and MCSH theory: a) gauge independent; b) gauge dependent. 

Fig. 6: Leading graphs contributing to the long distance behavior of < fifi* > in 

the broken phase of the AHM and MCSH theory: a) gauge dependent; b) gauge 

independent. 
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